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We show that one of exact solutions of the second-order Dirac-Coulomb equation were 
pseudo. In the corresponding original literature, it was considered that the mass of the 
electron with a Coulomb potential was position-dependent, but the obtained eigenvalues set 
was not the inevitable mathematical deduction of the given second-order Dirac equation, 
and the second-order Dirac equations were not the inevitable mathematical deduction of the 
given couplet first-order Dirac equation with the position-dependent mass of the electron. 
In the present paper, we obtain the correct solution of the introduced first-order differential 
equations. This new solution would be tenable only when the wave equation is correct, but 
there is not any experiment date to validate the so-called position-dependent of the electron 
in the Coulomb field. 
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I. INTRODUCTION 

An exact solution of the Dirac equation with the so-called position-dependent mass of the 
electron in the Coulomb field was introduced [1]. According to the article, in atomic units (H = 
tuq = 1), the spherically symmetric singular mass distribution was taken as follows 

m (r) = 1 + fi\ 2 It (1) 

where A is the Compton wavelength hjra^c = c _1 , and \i is a real scale parameter with inverse 
length dimension. The author presented some reasons why he introduced the position-dependent 
mass of the electron. We find the reasons were independent of any physical logic. For example, 
it was of that the rest mass of the particle (mo = 1) was obtained either as the asymptotic limit 
(r — ► oo), or the nonrelativistic limit (A — ► 0) of m (r), consequently, a possible interpretation for 
this singular mass term might be found in relativistic quantum field theory. It was even told of 
that it should also be noted that this position-dependent mass term has a relativistic origin as 
well since it was proportional to the Compton wavelength which vanishes as c — ► oo(equivalently, 
A -► 0). 
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Whereas we firstly query a question here, did any experiment ever show that the mass of 
electron in the Coulomb field is relative to the position? One can also find that many such 
suppositions in the corresponding published paper are very vexed. We are clear that there is not 
any consequence for the mo and relativistic result. Of course, at present, we have to transitorily 
avoid such questions and only check the corresponding mathematic deduction procedure. We show 
that the original solution of the second-order Dirac equation with the so-called position-dependent 
mass was incorrect, and the given second-order Dirac equation is not the necessary mathematical 
deduction of the given first-order Dirac equation with position-dependent mass of the electron in 
the Coulomb field. We introduce the correct exact solution of the original first-order differential 
equation only for further showing that the original solution includes many mathematical mistakes, 
and don't think the supposition of the position-dependent mass of the electron in the Coulomb 
field is correct. 



II. ESTABLISHED SOLUTION OF DIRAC EQUATION WITH 
POSITION-DEPENDENT MASS 

In order to solve the Dirac equation with the so-called position-dependent mass term in the 
Coulomb field, the spinor wavefunction was written as follows 

\[f(r)/r]<J-xim) 

where / and g are real radial functions, f is the radial unit vector, and the angular wavefunction 
with the spherical harmonic function Y t for the two-component spinor was written as 

= _J_ ( ^±m + l/2 Vr 112 

Xim V2TTT 1 Ty / LTrn + 1/2 Y r +1 ' 2 ' 



(4) 



then the following 2x2 matrix equation for the two radial spinor components was given 

l + A 2 ^-^ X(i-i) \ lair) 

+ -l + X 2 ^-e) \f(r) 

where e is the relativistic energy which is real. By introducing some signs and using the global 
unitary transformation, it was alleged that, the Schrodinger-like wave equation was finally obtained. 
The main steps are as follows. The Schrodinger-like requirement dictates that the parameter r\ 
satisfies the constraint 

Ch + Sk/X = ±Z (5) 
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where S = sin (Xrj), C = cos (A77) and —n/2 < A77 < vr/2, k = dbl, ±2, • • • . The solution of the 
constraint ([5]) gives two angles whose cosines are 



C=( fl 2 + K 2 /A 2 )- 1 



> 



(6) 



the equation Q is now transformed into the following 

C-e + (l±l)A 2 f A(-$ + ?-i) 
A(-S + ? + i) -C- £ + (lTl)A 2 f 




(7) 



where 7 = ^ -^/k 2 + A 2 (// 2 — Z 2 ) and 



cos 



A'/ 



sin 



An 



sin^ cos^ / V / 
equation ([7]) gives one spinor component in terms of the other as follows 

, A ( S 7 d 



whereas, the resulting Schrodinger-like wave equation becomes 

d 2 | 7 (7±1) | 2 Ze + y e 2 -l 



dr 2 



A 2 



4> T (r) = 



(8) 



(9) 



(10) 



the author comparing this equation with that of the well-known nonrelativistic Coulomb problem 

with constant mass and alleged that he had new discover for the relativistic spectrum 

2" 



1 + 



n + 1 + 1 



-1 



(n + L + iy 



1 + A 2 - 



Z 2 -n 2 
(n + i+iy 



(11) 



namely 



X 2 fiZ ± sj{n + 1 + l) 4 + A 2 (Z 2 - /i 2 ) (n + 1 + 1) 
{n + i + lf + {\Z) 2 



(12) 



where t stands for either one of the four possible alternative values in the following expression 
associated independently with cfy^ and ±k > 1. For <p + 



and for 



7 or 1 -+ -7 - 1, Z -► Ze + /x, £ -> (e 2 - l) /2A 2 



7 - 1 or t -> -7, Z -» Ze + /x, £ -> (e 2 - l) /2A 2 



(13) 



(14) 



It is similar to the above procedure, writing the paper about the relativistic quantum mechanics, 
many authors did not introduce their detail operation steps on how to obtain those necessary 
transition equation and their new mathematical result. They only alleged that those formulas in 
their paper are necessary deduction. We don't understand what the expression Q means? 
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III. EINGEVALUES-SET flU]) DISOBEY UNIQUENESS OF SOLUTION OF WAVE 

EQUATION 

For the same quantum system, it should have only one of eigenvalues set for any theory. The 
formula (|12p of the energy levels includes two eigenvalues set corresponding to different definitions. 
It is one of the mathematical contradictions of the articles. Consequently we cannot believe the 
formula (|12h is the real energy eigenvalues set in the Coulomb field. It seems that the author is 
not up on the method of finding the eigensolutions of the second-order differential equations with 
variable coefficients. Why don't we directly solve the second-order differential equation (jlOp now? 
Now, one can write the equation (jlOH in the separate form 



d 2 7(7+1) r, Ze+fi _ e 2 -} 

dr 2 r 2 r X 2 

d 2 , 7(7-1) n Ze+fj, _ e 2 -\ 

W 1 r 1 r * r "A 2 " 



<f>(r) = 
i/t(r) = 



(15) 



using the asymptotic solution <fi (r) ~ exp (—ry/l — e 2 j X^j , t/j (r) ~ exp (— r\/l — e 2 j X^j satisfying 
the boundary condition at r — > 0, we seek the formal solution 



\/i-4 \f 



.-a 



it easily obtained that 



(f)(r) = e a r M , tp (r) = e * r v (16) 



d4>(r) _ _ 

dr \ dr A 



lP~- e A 1 A — dF + ^^ u 



^flr ( dv Vl^l. 



dip(r) _ 

dr ° \ dr A 



(17) 



d^Kr) _ r -^^r ( d 2 v VT=^| d i_ £ 2 
~dP~ ~ e [dF 1 A dr~ + ^r V 



substituting (fl6|) and (fTT|) into (fT5|) . it educes that 



d 2 « _ ayi-eS du 7(7+1) ,, _ Z £lt +u _ n 

dF A dr « ^ r n — U 

^ _ 2^1-eg ^ _ 7(7-1) _ 9 Ze„+» _ n 
dF A dr P U Z r y — U 



(18) 



finding the power series solution of the above equations, it assumed that 

00 00 
u = Kr Su+n , v = J2 d n rSv+n (19) 



n=0 n=0 



so that 



t = E K + n) 6 n r s "+™- 1 

n=0 

oo 

& = £ («* + n) dnr^- 1 

n=0 

oo 

£ (s u + n)(s u + n- 1) 6 n r s 

n=0 
,9 00 

C£ — V ( a -L m\ ( g X«_1U T .s„+n-2 
n=0 



d 2 « 
dr 2 



r.Su+n — 2 



(20) 



X] 0„ + n) (s„ + n - 1) d n r 

ri=0 

substituting (fl~9|) and ([20j) into the equations (fT7|) . we have 

[{s u + n) (s u + n - 1) - 7 (7 + 1)] & n 



00 

n=0 I —2 



(s u + n - 1) + (Ze u + n) 



r 



s u +n— 2 



00 

£ < 

n=0 



[(s„ + n) (a„ + n - 1) - 7 (7 - 1)] d n 
^3 («„ + n - 1) + {Ze v + jj) 



(21) 



„s„+n-2 _ q 



finally we obtain the two recursive relation for the power series (Qj) 

[(flu + n) {s u + n - 1) - 7 (7 + 1)] - 2 
[(s„ + n) (s„ + n - 1) - 7 (7 - 1)] d n - 2 



(a« + n - 1) + (Ze u + fx) 



(a„ + n - 1) + (Ze,, + n) 



(22) 



using the initial value condition 6_i = 6-2 = • • • = 0, d-\ = d_2 = • • • r 
n = and substitute it into the above recursive relation, it educes that 

s u (s u ~ 1) - 7 (7 + 1) = 
s v (s v - 1) -7(7- !) = 



&n-l = 

d„_i = 
Oand 6 / 0, d ¥= 0, let 



(23) 



this gives 



Sui = 1 + 7, s u2 = -7; fl^l = 7, Si>2 = 1 - 7 



(24) 



in order to that the whole solutions ([2]) of the original equation satisfy the boundary condition, we 
have to choose 



hence 



u 



s u x = l + 7, s vi = 7 



n=0 



ra=0 



and 



2 OO 



1 ^ 00 



cf) (r) = e~^ Jtr ^ Kr 1+1+n , ijj (r) = e~^^ r ^ d n r^ + 



(25) 



(26) 



(27) 



n=0 



n=0 
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the power series must be cut off so that the whole wave function is limit at r — > oo. It is assumed 



that b nr ± 0, d nr / and b Hr+ i = b nr+2 = • • • = 0, d nr+ i 
recursive relation (|22l) . let n = n r + 1, we obtain 



(*« + n r ) + (Ze v + ^ 



n r +2 " 






0. According to the 



(28) 



it requires that 



^3 (i + 7 + nT ) + ( Z£u + /i) = 
( 7 + n r ) + (Ze„ + /i) = 
if we think little of, we would obtain the formal solution 



(29) 



£«2 



-A 2 Z Ai +^/(l+7+n r ) 4 +A2Z2[AV 2 + (l+7+™r) 2 ] 

A2z2+(l +7+ n r ) 2 

-A 2 Z/ i - v /(l+7+n r .) 4 +A 2 Z2[A 2 M 2 +(l+7+nr-) 2 ] 

A 2 Z2+(l+ 7 +n r ) 2 
-A 2 ZM+^/( 7 +n r ) 4 +A2Z2[ A 2^2 + ( 7+nr) 2j 

A2z2+( 7+ n r ) 2 
^A 2 Z At -y(7+n r ) 4 +A2Z 2 [A 2 ^ 2 + (7+n r ) 2 ] 
A 2 Z 2 +( 7 +n r ) 2 



(30) 



in form, these results as the inevitable deductions of the second-order differential equations (Q2 
include omnifarious logic problems. 

a) The solutions (|30l) are different from the formula (I12p , It shows that the original formula (|12() 
is incorrect for the second-order differential equations ([U). 

b) For the same quantum system described by the equations (JHJ), the four eigenvalues sets 
of the energy levels disobey the uniqueness of the solution of the differential equations. Which 
eigenvalues set is correct? 

c) Is the relativistic energy the positive number or the negative number? If we delete the 
negative energy solution, we also have two eigenvales set corresponding the positive energy. They 
also disobey the uniqueness of the eignesolutions set for the differential equation. 

d) It is the most serious that, the solutions (|30p implying the (|12p are the formal solution. 
Because of the definition 7 = k? + A 2 (fj, 2 — Z 2 ) given in the original article, only when k < 0, 



the quadratic equations with one unknown can have the solutions (j30l) . Whereas the k values 
constructed by Dirac are actually dbl, ±2, • • • , also given in the original articles. 

All of these problems are the mathematical and physical contradictions. Consequently, we don't 
think the original solution (|12p is not deceitful solution. 
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IV. SCHRODINGER-LIKE EQUATION (HOj IS BOGUS 

We even doubt that the second-order equation (fTUl) is the correct deduction of the original 
coupled first-order equation (j4j). By all appearances, from ((4|) to (jlOh . it actually introduces the 
transformation ([8]) with Q, namely 



C+e V A r ^ dr) V 

A C_S i 7 I d.^ 
C-e V A r dr ^ 




it gives the separate form 



C+elA rJ~drJZ C-e\ \ r dr )^ 



C osM +s i n M 
X _ C+e V A r*dr) v C — e V A r rfr / 



cos^i+sin^i 



(31) 



(32) 



(33) 



these transformations cannot translate the equations (JH) namely 

(l + A 2 ^-e) 5 (r) + A(f- J:)/(r) = 
Mf + if) 9(0 - (-1 + A 2 ^ - e) / (r) = 

into any Schrddinger-like second-order equations, and the second-order equation (|1U|) give in the 
original paper is not correct. 



V. CORRECT SOLUTION OF THE CORRESPONDING FIRST-ORDER EQUATION 

Because there has been not any experiment data to approve the imagination, we don't think 
that the mass of electron in the Coulomb field is really dependent on position. Consequently, in 
principle, it should be meaningles to find the correct solution of the equation @ or (|33p . When only 
looking from a mathematical point of view, we give the correct exact solution of the equations (J1J) or 
(I33[) . It is well known that directly solving the original coupled first-order equation @ namely (|33|) 
is simpler out and away than translating it into the so-called Schrodinger-like equation to obtain 
the exact solution. One firstly note the behavior of / (r) and g (r) for g —* oo, since neglecting the 
terms proportional to 1/r the differential equations ([33]) read 

(1- £ ) 5 -a|J = 0, A^ + (l+e)/ = (34) 



it follows immediately that 



d 2 f _ d?gl-e* 

d^ + — f ~°> d^ + — 9 ~° (35) 
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there will be different solution with the different fields of definitions of the e. However, it is 
considered all along that the relativistic energy in the Coulomb field satisfy the condition < E < 
tooc 2 namely E < 1. According to ()13f) and (j!4[) . it only gives 1 < e 2 < 2A 2 + 1. For the moment, 
it is considered that e 2 > 1. We obtain the asymptotic solutions of the equations (f33j) 



and the exact solution of the equations 



, g ~ e 
take form 



(36) 



they give 





v 7 ^" 


-i 


f = e~ 


A 


T 

u 


&-e~ 


A 


—r du 


dr ~ e 




dr 


d 9 _ „-- 




—r dv 




X 


dr 
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du Ve 2 — 1 



(37) 



A U 



dv _ x/e^T 



(38) 



substituting for equations (|33p . we have 



A f - if + V^T)u + 



\dv , 
A dF + V r 



Xk 



v + 



(s - 1) 



A 2 (Z+ M ) 



v = 
it = 



(39) 



„cr+n 



finding the power series solution, put 

oo oo 
n=0 n=0 

substitute it into the above equations, we obtain 

OO / \ OO r 2 / v \ ~ 00 

A E (ff + n) dnr^- 1 - ( f + ) £ d n r° +n + (e - 1) - £ ^r ff+n 



(40) 



n=0 

oo 



n=0 

oo 



A £ (a + n) 6 n r CT +"- 1 + ( M - E b n r° +n + (e + 1) - 

n=0 ^ ' n=0 L 

it predigests that 

OO i- 

E A (a + n - k) d n - v/^Tdn-i + (e - 1) 6 n -i - A 2 (Z + //) 6 r , 

n=0 L 

OO r 



n=0 

oo 

E 4r ff+n 

n=0 



(41) 



E A 2 [Z - n) d n -(e + 1) d n -i - X(a + n + k)b n + \/e 2 - 16 n -i 

71=0 L 

so the coefficients of the power series satisfy the recursive relations 



„<T+n— 1 



„<7+rt— 1 



(42) 



A (a + n - fc) d„ - \/e 2 - ldn-i + (e - 1) &n-i - A 2 (Z + p) b n = 
A 2 (Z -(j)dn-(e+ 1) d n -i - A {a + n + k) b n + Ve 2 - 16„-i = 



(43) 
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Solving the above recursive relations and using the initial conditions that bo 7^ 0, 6_i = 6-2 
= and do = 0, d-\ = = • • • = 0, put n = in (H3|) . we obtain 



A (ex - A;) d - A 2 (Z + n) b = 
A 2 (Z - n) d - A (0- + k) b = 

it requests that the determinant of coefficient is equivalent to zero, so that 



(44) 



A (a - k) -A 2 (Z + n) 
X 2 (Z-fl) -\(a + k) 







(45) 



it educes that value of the index a = ±\Jk 2 + A 2 (Z 2 — fi 2 ), since the wave function must be limited 
at r — > 0, it can be only taken 



a = yjk 2 + A 2 (Z 2 - ll 2 ) 



(46) 



on the other hand, the formal whole wave function should be combined by expressions (|2|), (|3|), 
(J37D, & and flUD, that is to say 

/ „-„-^=ir V h r V k2 + x2 ( z2 -^ 2 )+n-lJ \ 
Z_i Aim 



n=0 

E 2_! OO 



V 



(47) 



n=0 



the boundary condition at r — ► 00 requests that the parts of the power series in the formal whole 
wave function must be cut off. It is assumed that b rir ^ 0, d Ur ^ and 6 nr +i = &n r +2 = • • • = 0, 
dnr+i = dn r +2 = • • • = 0, according to the recursive relations (fl3]) . put n = n r + 1, we obtain 



(48) 



-V^ld nr + (e-l)b nr = 

- (£ + 1) dn r + \/^T6n r = 

it indicates that the e should take some special values so that b nr+ \ = and dn r +i = 0. Note that 
the recursive relations (|43[) . Multiplying the first relation by \Je 2 — \ and multiplying the second 
relation by e — 1, we have 



A (a + n r — k) ye 



I d r . 



1) d 



n r — 1 



-A 2 (Z + fi) V^lbn, + (e - 1) V?^!^-! = 



A 2 (Z - /x) (e - 1) dr. 



1) d nr . 



(49) 



-A (a + n r + k) (e - 1) 6 nr + (e - 1) v 7 ? 



16^ 



this deduces that 



(a + n r + fc) (e - 1) - A (Z + /x) Ve^T^ 
(a + n r - k) y/e 2 - 1 - A (Z - fi) {e - 1) 



(50) 
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combing this relations and one of the system of the recursive relations (|48p , we obtain 



(cr + n r ) \J e 2 - 1 = AZe + Xfi 



(51) 



it gives that 



£ = 



A 2 



//Z ± (<r + n r ) \j {(J + n T ) 2 + A 2 (/U 2 -Z 2 ) 



(52) 



(<j + n r ) 2 - A 2 Z 2 



we find that it is different from not only the formula (jlip given in the original articles but also the 
correct formula (|30p of the second-order differential equations. This shows that the original exact 
solution of the Dirac equation with position-dependent mass is the pseudo solution. 



We have used the basal knowledge of the differential equation with variable coefficient to show 
that, almost every step of the mathematical operation in the mentioned paper for introducing the 
Dirac theory with the position-dependent mass of the electron in the Coulomb field is incorrect, 
and the energy eigenvalues so the eigen-wave-function given in the original paper are pseudo. 
Only in a mathematical signification, we give the correct solution and eigenvalue set of the first- 
order differential equations. However we don't regard the formula (I52p of the energy-levels as the 
necessary result of the development of quantum mechanics. Because the position-dependent mass 
of the electron in the Coulomb field make us discredit its authenticity. 

Some other papers also alleged that they found the exact solution of the Dirac equation for a 
particle with position-dependent mass [2], which might be useful in the study of the corresponding 
non-relativistic problem as a reference result. It was even considered that the next terms, which 
they have neglected in this work (in particular the dipolar one) and which are responsible for the 
super-fine structure of the energy spectrum, can be taken into account by means of standard per- 
turbation theory. Although the corresponding paper constructed the second-order Dirac equation 
by using unconventionality methods which also can be seen in the some papers published 25 years 
agoQ]. However, is this theory correct? Without solving the corresponding differential equation, 
one can find some methods to directly conclude that some other papers about the Dirac theory are 
not correct. 
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